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Recently, the chiral-induced spin selectivity in molecular systems has attracted extensive interest 
among the scientific communities. Here, we investigate the effect of the gate voltage on spin-selective 
electron transport through the a-helical peptide/protein molecule contacted by two nonmagnetic 
electrodes. Based on an effective model Hamiltonian and the Landauer-Buttiker formula, we calcu¬ 
late the conductance and the spin polarization under an external electric field which is perpendicular 
to the helix axis of the a-helical peptide/protein molecule. Our results indicate that both the mag¬ 
nitude and the direction of the gate field have a significant effect on the conductance and the spin 
polarization. The spin filtration efficiency can be improved by properly tuning the gate voltage, 
especially in the case of strong dephasing regime. And the spin polarization increases monotonically 
with the molecular length without the gate voltage, which is consistent with the recent experiment, 
and presents oscillating behavior in the presence of the gate voltage. In addition, the spin selectiv¬ 
ity is robust against the dephasing, the on-site energy disorder, and the space angle disorder under 
the gate voltage. Our results could motivate further experimental and theoretical works on the 
chiral-based spin selectivity in molecular systems. 

PACS numbers: 87.14.E-, 87.15.A-, 85.75.-d, 73.63.-b 


I. INTRODUCTION 

Molecular spintronics, which manipulates the elec¬ 
tron spin transport through organic molecules, has 
been receiving more and more attentions, because of 
the potential applications in storage and processing of 
informationi 1_ ' A growing number of pioneering re¬ 
searches have established organic materials as a vi¬ 
able platform for spin-filter devices. The experimen¬ 
tal evidence of room-temperature spin-polarized injec¬ 
tion and transport through prototypical organic semicon¬ 
ductors was demonstrated for the first time£ An orig¬ 
inal spin-valve device was designed in which a single- 
walled carbon nanotube is laterally coupled to single¬ 
molecule magnets through supramolecular interactions^ 
Spin-selective carrier transmission was demonstrated in 
a nonmagnetic system which is composed of carbon 
nanotube-DNA hybrid j 10 : 11 Both self-assembled mono- 
layers of double-stranded DNA (dsDNA) deposited on 
gold substrate^ and single dsDNA molecules connected 
by two electrodes^ presented high spin polarization at 
room temperature. However, spin-polarized effect was 
not observed in single-stranded DNA monolayers— Later 
on, spin-selective effect of electron transmission along 
bacteriorhodopsin—an a-helical protein—embedded in 
purple membrane physisorbed on gold and aluminum 
substrates was demonstrated— Some works theoretically 
investigated the spin transport properties of the chiral 
molecular systems — ~-° A model Hamiltonian, includ¬ 
ing the small environment-induced dephasing, the weak 
spin-orbit coupling (SOC), and the helical symmetry, 
was proposed to perfectly rationalize the quantum spin 
transport through the dsDNA and the single-stranded 
DNA , 2 - i 22 as well as the a-helical protein, and meanwhile 


explained the contradictory results between the protein 
and the single-stranded DNA—. On the basis of self- 
assembled monolayer of a-helical polyalanine adsorbed 
on gold, a device was presented and indicated the abil¬ 
ity to produce spin-based device without a permanent 
magnet— In addition, it was demonstrated that the elec¬ 
trons transmitted through Photosystem I which is mainly 
composed of protein complex were highly spin-selective 
and the spin polarization was temperature-dependent— 
Very recently, the light-controlled ability on spin filtra¬ 
tion through the bacteriorhodopsin D96N mutant was 
observed— The spin selectivity of electron was measured 
through the monolayers of oligopeptides and increased 
with the increase of the molecular length in the investi¬ 
gated length rangei 27 

Recent experimental and theoretical researches on 
spin-dependent electron transport through the DNA and 
protein molecules have given rise to a prominent im¬ 
provement of molecular spintronics. Electron transport 
through the dsDNA molecules presented high spin polar¬ 
ization. Meanwhile, the dsDNA molecules could act as a 
field-effect transistor under a gate voltage.— Then a the¬ 
oretical investigation on the influence of gate voltage on 
spin transmission along the dsDNA molecules was per¬ 
formed, revealing that spin polarization showed strong 
dependence on the magnitude as well as the direction 
of gate voltage and could be significantly enhanced by 
tuning the gate voltage— As for the a-helical protein, 
the transmitted electrons also exhibited the ability of 
spin filtration. In the presence of gate voltage, protein 
field-effect transistors were also reported and a model 
for transport was proposed. 30,31 Based on the similarly 
unique helical structure, spin-selective properties of the 
a-helical protein were compared with that of the DNA 
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molecules and then a question came to us that whether 
the gate voltage has also such a intense effect on the spin 
transport along a-helical protein, just like along DNA 
molecules. In addition, we also wanted to investigate 
how the on-site energy disorder and space angle disorder 
affect spin polarization of the peptide in the presence of 
gate voltage. 

Herein, we report on a method to regulate the spin- 
dependent electron transport along the a-helical protein 
molecule connected by two non-magnetic electrodes in 
the presence of gate voltage, which gives rise to an ex¬ 
ternal electric field perpendicular to the helix axis of the 
molecule, as illustrated in Fig. [T| On the basis of an 
effective model Hamiltonian and the Landauer-Biittiker 
formula, the conductance and the spin polarization are 
calculated. Our results indicate that the spin filtration 
efficiency shows strong dependence on the magnitude as 
well as the direction of the gate voltage. One can im¬ 
prove the spin polarization by properly tuning the gate 
voltage, especially in the case of strong dephasing regime, 
such as the high experimental temperature. Both the 
conductance and the spin polarization versus the protein 
length show oscillating behavior when the gate voltage 
is employed. We also find that the spin polarization of 
the peptide is robust against the dephasing, the on-site 
energy disorder and the space angle disorder under the 
gate voltage. Therefore the a-hclical protein-based de¬ 
vice could be a more efficient spin filter by properly tun¬ 
ing the gate voltage. 

The rest of the paper is organized as follows. In Sec. El 
the calculation model and the method are presented. In 
Sec. mu the effect of the gate voltage on the conductance 
and the spin polarization is shown. And then, we investi¬ 
gate the influence of the on-site energy disorder and the 
space angle disorder. Finally, the results are summarized 
in Sec. EU 

II. MODEL AND METHOD 

The spin transport properties of the a-helical protein 
can be simulated by the Hamiltonian^ 
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FIG. 1: (color online), (a) Schematic view of the a-helical 
protein under external electric field E g which is perpendicular 
to the helix axis (z-axis). The parameters are radius R = 0.25 
nm, Euclidean distance h = 0.41 nm, space angle 9 1 = 0.37 
rad, twist angle A<p = 5n/9, and stacking distance Ah = 0.15 
nm between two nearest neighbour amino acids. The amino 
acids are labelled by Arabic numerals from the bottom up. 
(b) Projection of the five bottom amino acids and the electric 
field into the x-y plane. Here, (fio is the angle between the 
direction of the electric field and the positive x axis direction. 


The first two terms are the Hamiltonian of the a- 
helical protein whose length is N, with cj) = (c^, c^) 
and c n = (c n -|-, c n 4,) T being the creation and annihi¬ 
lation operators, respectively. e n is the on-site en¬ 
ergy, t n j = j s thg hopping integral be¬ 

tween two neighboring sites n and n + j, p n j = 
SiCOs(ifi~j)e^ ( - lnj ~ ll ^ lc is the corresponding SOC, and 
<r n j = (&x sin <£+■ - a y cos (finj ) sin 9 nj + a z cos 9 nj 22- 

Here, l nj = y / [2i?sin(( / 5“ J )] 2 + ( h n+j - h n ) 2 is the Eu¬ 
clidean distance between sites n and n+j , = ( ip n +j ± 

<p„)/2, and 9 n j = arccos[2A sin(<p“j)/?„.,■] is the space an¬ 
gle, where the parameters ( R , ip n , h n ) are the three cylin¬ 
drical coordinates of site n. (J x ^ y ^ z are the Pauli matrices, 
l c is the decay exponent, and si is the renormalized SOC. 
In the absence of any disorder, l n j for j = 1 is reduced 
to l± = -\/[2I?sin(A<^/2)] 2 + (Ah) 2 and correspondingly 
t n j for j = 1 is reduced to t\ which denotes the hopping 
integral between the nearest neighbor (NN) sites, with 
Ap and Ah being the twist angle and the stacking dis¬ 
tance between the NN sites, respectively. The third and 
the fourth terms are the Hamiltonians of the left and 
right semi-infinite real electrodes, respectively. The next 
one represents the couplings between the a-helical pro¬ 
tein and the two real electrodes. Finally, the last term 
(1) denotes the Buttiker’s virtual electrodes and their cou¬ 
plings to each site of the molecule, which is introduced to 
describe the dephasing processes caused by the electrons’ 
inelastic scatterings with the electrons, the phonons, the 
counterions, and the adsorbed impurities j 32 ' 33 Actually, 
previous works have clearly indicated the presence of the 
phase-breaking processes in the proteins^ - — 

When the a-helical protein is subjected to an external 
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electric field which is perpendicular to its helix axis (see 
Fig. ID, the on-site energy at the nth site will be changed 
into the following form: 

£ n = + eE g Rcos(ip n + 4> 0 ), (2) 

where is the on-site energy without the external elec¬ 
tric field and e is the elementary charge. E g is the per¬ 
pendicular external electric field and thus the gate volt¬ 
age across the a-helical protein molecule is 2V g = 2E g R. 
The phase <j> o, which is the angle from the direction of 
the external electric field to the positive direction of x- 
axis, shows the orientation of the gate voltage relative to 
the helical molecule, as seen in Fig. |T| In order to adjust 
(j) o, the helical molecule could be rotated with the direc¬ 
tion of its helix axis fixed. Eq. © shows that the gate 
voltage harmonically regulates the on-site energies along 
the helical strand and introduces periodic change of each 
site, due to the intrinsic helical structure of the protein. 
Undoubtedly, such adjustment will affect the electronic 
structure of the a-helical molecule, which could make a 
significant effect on both the conductance and the spin 
polarization. The magnitude of the gate voltage chosen 
in this paper is the order of 0.1 V, where the external 
electric field is much smaller than the internal one gen¬ 
erated by the nuclei of the protein molecules and hence 
its effect on SOC may be negligible. 

From the Landauer-Biittiker formula, the current in 
the c/th real or virtual electrode with spin s 4- can 
be described as I qs = (e 2 /h)J2 m , s ' T q s,ms'(V m - V q ), 
where V q is the voltage of the c/th electrode and T qS)Tns ' 
is the transmission coefficient from the mth electrode 
with spin s' to the c/th electrode with spin s. With the 
boundary condition that the net current across each vir¬ 
tual electrode is zero, the voltage V q of the virtual elec¬ 
trodes can be derived by applying a small bias between 
the real electrodes with Vl = 14 and Vr = 0. Then 
the conductances of the right real electrode for spin-up 
(Gy) and spin-down (Gy) electrons can be written as 
G s = (e 2 /h) X)m,s' TRs,ms' V m /Vb. Finally, the spin po¬ 
larization is defined as: 


G t - G1 
~~ Gy ± Gy ’ 


and the averaged spin polarization is: 

m = hL p - dE - 


( 3 ) 

( 4 ) 


Here, fl denotes the lower energy band of E < E c and E c 
is the “band center”, where the number of the electronic 
states below E c is equal to that above E c . 

For the a-helical peptide, the structural parameters 
are the radius R = 0.25 nm, the twist angle A ip = 5 tt/9, 
and the stacking distance Ah = 0.15 nm. We take the 
NN hopping integral t\ as energy unit. The values of 
aforementioned parameters are chosen as, the molecu¬ 
lar length N=30, the on-site energy = 0 without 
loss of universality, the renormalized SOC parameter 


Si = 0.12ti, and the decay exponent l c = 0.09 nniS 
For the real electrodes, the retarded self-energy can be 
numerically derived from t m = 4ti and r = 2t±JH For 
the virtual electrodes, the dephasing strength is set to 
= 0.02fi. The values of all above-mentioned param¬ 
eters will be used throughout this paper except for spe¬ 
cific annotation. We also investigate the spin transport 
through the a-helical protein under the on-site energy 
disorder and the space angle disorder, as illustrated in 
Fig. [3 The significant effect of the gate voltage on spin 
transport along the a-helical protein molecule could be 
observed in a wide range of model parameters. 


III. RESULTS AND DISCUSSIONS 

As a comparison, we first consider the spin transport 
through the a-helical peptide with length N = 30 in 
the absence of the gate voltage. Figure [3(a) shows the 
spin-up conductance Gy (red-dashed line), the spin-down 
one Gy (green-dashed line), and the spin polarization P s 
(black-solid line) with V g = 0 and = 0.02ti. Although 
the dephasing strength is smaller than previous work, 23 
similar results can also be observed in the energy spec¬ 
trum. For instance, there exist several sharp peaks in 
the curves of Gy-U and G'y-U; the “band center” E c is 
shifted toward lower energy, i.e., E c < 0; except for the 
band center E c at which P s = 0, the spin polarization is 
nonzero and can achieve ±26.7%, which is in accordance 
with the experimental resultj 14 ' 27 



FIG. 2: (color online), (a) Energy-dependent spin-up conduc¬ 
tance Gy (red dashed line), spin-down one Gy (green dashed 
line), and spin polarization P s (black solid line) for the a- 
helical peptide without the gate voltage, i.e., V g = 0. (b) Gy 
and (c) P s for the a-helical peptide with different values of 
V g by fixing <j >0 = 0.47T. The case of V g = 0 is also displayed 
as a comparison. The legend in (c) is for both (b) and (c). 
(d) P s for several values of with V g = 0.8. 
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We then investigate the spin transport properties of the 
a-helical peptide under the gate voltage. Figures Ob) 
and He) display, respectively, the spin-up conductance 
G-j- and the spin polarization P s with several values of 
V g by fixing </>o = 0.47T, while the other parameters are 
the same as Fig. 0a). Here, the different colors of curves 
denote different V^’s. In the presence of the gate voltage, 
the period of the on-site energies is changed to be eigh¬ 
teen amino acids, leading to the appearance of multiple 
potential barriers and wells within each period. As a re¬ 
sult, the conductance is declined by increasing V g , and 
several deep transmission valleys could emerge around E c 
and will develop into gaps in the case of larger V g [see the 
blue- and green-dashed lines in Fig. [2jb)] . Besides, the 
transmission spectrum can move toward both lower and 
higher energies by increasing V g , since the range of the 
on-site energies (—eV g , eV g ) increases with V g . Although 
the transmission ability of the peptide is weakened by 
V g , its spin filtration efficiency is robust against V g and 
can even be enhanced by V g [Fig. 0c)]. For instance, 
the maximum of P s is 26.7%, 28.1%, 29.9%, and 20.2%, 
respectively, by increasing V g from 0, 0.4, 0.8, to 1.2. 
Meanwhile, for E < E c , the position of the peak becomes 
closer to E c \ for E > E c , the position of the valley (peak) 
remains the same as that of V g = 0. This phenomenon is 
different from the dsDNA molecules^ What’s more, it is 
interesting that the sign of the spin polarization P s can be 
reversed in the case of large V g [see the green-dashed line 
in Fig. [2jc)] . Figure 0d) shows P s for different values of 
<j >o by fixing V g = 0.8. Although the positions of the peak 
and the valley almost keep still by changing <f> o, the spin 
filtration efficiency of the peptide considerably depends 
on </>o and can be improved by properly modulating 4>o. 
For instance, when E = — 0.44fi, P s is 29.9%, 26.1%, 
20.0%, and 30.0%, respectively, by varying </> 0 from 0.47T, 
0.87T, 1.2-7T, tO 1.67T. 

Figures 0a) and 0c) plot, respectively, the conduc¬ 
tance G -(- vs the phase <f>o and the spin polarization P s 
vs the phase </>o with four values of the gate voltage V g . 
It is clear that in the absence of the gate voltage, both 
G -f and P s are independent of 0 O . When the gate voltage 
is applied, both the curves of G^-(f >o and P s -4 >o present 
oscillating behaviors. By inspecting Fig. 0a), it can be 
seen that for relatively small V g (V g = 0.4,0.8), the curve 
G^-(f>o displays significant oscillation phenomenon with 
three distinct peaks; while for large V g (V g = 1.2), the 
curve Gf-4>o possesses several peaks with quite small os¬ 
cillating amplitude, because of the strong gating effect 
at large V g . As compared with Gf-<j> o, the oscillating 
amplitude of the curve P s ~4 >o is always very big for all in¬ 
vestigated values of V g and the spin polarization is quite 
large (Fig. 0c)). One can see that the spin filtration ef¬ 
ficiency of the a-helical peptide at V g = 0.4 and 0.8 is 
always larger than that without the gate voltage. Even 
in the case of large V g , there still exist several intervals 
of </>o at which P s is larger than that of V g = 0. There¬ 
fore, the spin filtration efficiency of the a-helical peptide 
could be drastically enhanced by properly adjusting the 



FIG. 3: (color online), (a) Gf vs <j >o and (c) P s vs 4 >o for the 
peptide with different values of V g . (b) Gf vs V g and (d) P s vs 
V g for the peptide with several values of </>o- Two-dimensional 
plot of (e) P s and (f) (Ps) vs V g and (f>o- It is obvious that P s 
and (P 3 ) have strong dependence on V g as well as <j> o. Here, 
the Fermi energy is E = — 0.41ti. 


direction of the gate field. 

We then show the conductance G -j- and the spin polar¬ 
ization P s as a function of the gate voltage V g with four 
different values of </>o, as illustrated in Figs. [3)(b) and (d). 
One can note that the dependence of both G -j- and P s 
on V g is more concise, which is different from the curves 
Gf-(j)o and P s -(f>o- The conductance G-j- is almost inde¬ 
pendent of V g for V g < 0.3 and decreases sharply with 
increasing V g for 0.3 < V g < 1.1 (Fig. 0(b)). And G-j- 
tends to zero for V g > 1.1, since the transmission valley 
could appear around the energy E = — 0.41fi due to the 
strong gating effect. In contrast, the dependence of P s 
on V g is not monotonic and there exists a turning point 
Vg in the curve P s ~V g . P s increases with V g at first for 
V g < V g and is gradually declined by further increasing 
V g . This phenomenon is irrespective of <f> 0 , although the 
turning point may depend upon the direction of the gate 
field. Furthermore, P s of the a-helical peptide is rela¬ 
tively large for very large gate voltage, e.g., V g > 1.5. 
This indicates that the spin selectivity of the peptide 
could be enhanced by changing the gate voltage and is 
robust against the strong gating effect. 

In the following, we further investigate the spin polar- 
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FIG. 4: (color online), (a) P s vs length N with Yd = 0.02ti 
and (b) P„ vs N with Yd = 0.06ti for the peptide/protein at 
E = — 0.41ti. (c) Averaged conductance (G-f) vs N and (d) 
(P a ) vs N for the peptide/protein with Fd = 0.02fi. Here, 
0o = 0.87T, the maximal length N = 80, and the curves of 
different colors represent different V g . 


ization P s of the ct-helical peptide with the magnitude 
and the direction of the gate field in a wider parameter’s 
range, as illustrated in Fig. EHe). One can see that P s 
increases with V g at first and then is declined by further 
increasing V g , for almost all values of 0o . And there exist 
three extensive domains at which P s exceeds 20% and is 
much bigger than that of V g =0. The approximate range 
of these three domains is [0.7,1.1] for V g and [0.37T, 0.57t], 
[0.97t,7t], [1.57T, 1.87t] for 0 O . Besides, one can identify 
some other features. (I) For V g < 0.3, P s is nearly inde¬ 
pendent of 0o- (II) For 0.3 < V g < 1.1, the dependence 
of P s on 0o exhibits three peaks. (Ill) For V g > 1.1, the 
behavior of P s versus 0o becomes complicated and has 
multiple turning points. The above results reveal that the 
spin polarization of the peptide is strongly dependent of 
the magnitude as well as the direction of the gate field, 
and could be significantly increased by properly tuning 
the gate voltage. Figure E]T) displays the averaged spin 
polarization (P s ) of the peptide as functions of V g and 0o- 
It can be seen that (P s ) monotonically declines with the 
increase of V g , irrespective of 0o- {Ps) presents oscillat¬ 
ing dependence on 0o by fixing V g . With the increase of 
V g , the peak width of (P s )-0o decreases and some small 
peaks appear. 

Figures Oja) andQJb) show the length-dependent spin 
polarization for four different values of V g with = 
0.02ti and = 0.06ti, respectively. One can see the 
following features which are similar for different values 
of r d . (1) In the absence of the gate voltage, P s in¬ 
creases monotonically with N £2- For small the de¬ 
pendence of P s on N is almost linear and the rising slope 
decreases very slowly with increasing N, which is in ex¬ 
cellent agreement with recent experiments 27 While for 


relatively large P s increases very fast with N for 
short molecular length and the increasing speed slows 
down for long molecular length, since larger dephasing 
strength leads to faster memory loss of the electrons. (2) 
When the gate voltage is implemented, the curve P s - 
N displays oscillating behavior due to the gating effect. 
One could note that the value of the same figurate peak 
increases quickly with N at first and afterwards is sup¬ 
pressed with the further increase of N, which is more 
obvious for large V g . (3) The oscillating period is length 
18, which corresponds to the same figurate peak value 
and derives from the fact that the twist angle of the a- 
hclical protein is A ip = 57r/9. However, the oscillation 
behavior of the peptide/protein is much more compli¬ 
cated and several peaks could appear within a period, 
because a period is composed of five helical circles. This 
is different from the dsDNA molecule,— where only one 
peak exists in a period, since a period is composed of 
one helical circle in the dsDNA molecule. (4) The oscil¬ 
lating amplitudes of the P s -N curves for V g = 0.8 and 
V g = 1.2 are much larger than that for V g = 0.4. (5) 
For all investigated values of the gate voltage, there al¬ 
ways exist specific length regions, where P s under the 
gate voltage is larger than that without the gate voltage. 
The range of these specific length regions decreases with 
increasing V g . For instance, P s at V g = 0.4 is larger than 
that at V g = 0 for almost all the length region. P s at 
V g = 1.2 is larger than that at V g = 0 only for short 
molecular length. These results further indicate that the 
spin selectivity of the peptide/protein could be improved 
by introducing an external electric field in a quite wide 
length range, especially for short peptide. 

Figures 0][c) and[4][d) plot the averaged spin-up con¬ 
ductance (G-j-) versus N and the averaged spin polar¬ 
ization (P s ) versus N, respectively, with = 0.02ti 
for several values of V g . In comparison with the case 
of V g = 0 that both (G-f-) and (P s ) vary smoothly with 
N, the behaviors of {Gf)-N and (P s )-N under the gate 
voltage are oscillating and the oscillating amplitude is al- 



FIG. 5: (color online). Length-dependent P a for the pep¬ 
tide/protein with three values of 0o at E = — 0.41ii and 
Vg = 0.8. 
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most the same for different values of V g and N. It can be 
seen that the peak value of {G^)-N declines fast with N 
at first and is suppressed for long molecular length, espe¬ 
cially in the case of large gate voltage. Contrary to the 
dependence of (Gf) on N, (P s ) increases monotonically 
with N in the absence of gate voltage, which is consistent 
with the experiment i2Z While for V g 7 ^ 0 in the ( P s )-N 
curve, the value of the same figure peak increases with N 
at first and then tends to saturation by further increasing 
N, especially in the case of large V g . The larger V g cor¬ 
responds to the smaller values of (G-j-) and (P s ), which is 
independent of N, because of the gating effect. What’s 
more, both the averaged conductance (G-j-) and the aver¬ 
aged spin polarization (P s ) are still comparatively large 
for N = 80 and V g = 1.2. 

Figure [5] displays the length-dependent spin filtration 
efficiency of the a-helical peptide/protein for different 
values of the phase <f >0 at = 0.02ii. P s oscillates 
by increasing N in the presence of the gate voltage, for 
whatever the value of fa ■ The curve P s -N for various fa 
presents two similar oscillating features as that of the P s - 
N curve for different V g . Firstly, the same figurate peak 
shows that the period of the length is 18 amino acids. 
Secondly, the value of the same figurate peak increases 
with N at first and then is suppressed by further increas¬ 
ing N. Although the oscillating amplitude of the curve 
P s -N nearly shows independence on <j> 0 , the position of 
the same figurate peak moves towards larger N with fa 
at first and then is shifted towards smaller N by the fur¬ 
ther increase of <j> 0 . For instance, one local maximum of 
the corresponding peak locates at N = 68 , 74, and 73, 
by increasing fa from 0 . 67 T, 1.27r, to 1 . 87 T. It can also be 
seen that there exists a quite wide length range where 
P s varies drastically with fa by fixing N. The results 
imply that the spin polarization of the peptide/protein 
for various length is (^-dependent and could be tuned by 
adjusting the direction of the gate field. 

The dephasing occurs inevitably in the experiment and 
its strength may rely on various experimental conditions, 
such as the temperature. Here, we consider the influence 
of the dephasing on the spin transport through the a- 
hclical peptide under the gate voltage. Figures [(Ja) and 
[ 6 tT) plot the spin polarization P s and ( P s ), respectively, 
as a function of with the strength up to = t\ for 
four values of V g . It is clear that there exists a crossover 
of in each curve of P s -Td and (P s )-Td, although the 
specific value of this crossover may depend on the gate 
voltage. Both P s and (P s ) increase quickly with in 
the weak dephasing regime and then are slowly declined 
by further increasing Td in the strong dephasing regime. 
This phenomenon does not depend upon the value of the 
gate voltage and originates from the two competitive ef¬ 
fects of the dephasingi^ On the one hand, the dephasing 
promotes the openness of the system because of the cou¬ 
pling to the Biittiker’s virtual lead and produces the spin 
asymmetry. On the other hand, the dephasing gives rise 
to the memory loss of the electrons and shrinks the spin 
polarization. When the former effect dominates, the spin 
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FIG. 6: (color online), (a) P s vs dephasing strength Vd at 
E = —0.41ti and (b) ( P s ) vs I’d for the peptide with various 
values of V,, by fixing tj> 0 = 0.87T. 


polarization could be enhanced by increasing when 
the latter effect prevails, the spin polarization will be de¬ 
clined. Besides, one can see from the curve (P s )-Td that 
the position of the crossover is shifted toward larger T^ by 
increasing V g and (P s ) at the crossover is declined. How¬ 
ever, in the strong dephasing regime where the spin polar¬ 
ization decreases with increasing the absolute value 
of the decreasing slope is declined by increasing V g . This 
indicates that the spin filtration efficiency could be more 
robust against the dephasing in the region of large V g . 
In other words, both P s and (P s ) could be enhanced by 
the gate voltage in the strong dephasing regime. There¬ 
fore, the gate voltage could enhance the spin filtration 
efficiency of the a-helical peptide in a wide dephasing 
range, especially in the case of the strong dephasing. 

Distinct amino acids that form the a-helical pep¬ 
tide/protein could lead to the on-site energy disorder. 
Next, we study the influence of the on-site energy disor¬ 
der on the spin transport through the peptide molecule 
in the presence of the gate voltage. A stochastic variable 
w n is added to each e n to simulate random on-site en¬ 
ergy disorder with w n uniformly distributed in the region 
[-W/ 2, IT/2] and IF is the disorder strength. Figurc[7]]a) 
plots the spin polarization P s versus the energy E for 
five different values of IF with V g = 0.4 and fa = 0.87T. 
One notes that the positions of both the peak and the 
valley move towards higher energy with increasing IF 
and the corresponding P s is usually declined, due to the 
disorder-induced Anderson localization effect. However, 
P s at some energy regions could be enhanced by increas¬ 
ing IF and is robust against the on-site energy disorder 
in the presence of the gate voltage. This statement is fur¬ 
ther demonstrated in Fig. [3b), where (P s ) is displayed 
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FIG. 7: (color online), (a) P s vs E under the on-site energy 
disorder with several strength W and (c) P„ vs E under the 
space angle disorder with several strength D for the peptide 
by fixing V g = 0.4 and <j >o = 0.87T. (b) (P a ) vs on-site en¬ 
ergy disorder strength W and (d) (P s ) vs space angle disor¬ 
der strength D for the peptide with different V g at </>o = 0.87T. 
All the data are performed for a single disorder configuration 
and similar results can be also obtained with other disorder 
configurations. 

as a function of W for several values of V g . It clearly 
appears that (P s ) could be increased by increasing W in 
the weak disorder region and is still considerably large in 
the strong disorder region. 

Each amino acid will fluctuate around its equilib¬ 
rium position, which can be simulated by introducing 
a random variable d n to each 9 n \ between the nth and 
(n + l)th sites. d n is evenly distributed within the range 
[— D/2, D/2] with D the disorder strength of the space 
angle. Herein, we assume that all of the amino acids 
locate at the side surface of the cylinder whose radius 
R is constant and the Euclidean distance l n \ is always 


the same as Then, the stacking distance and the 
twist angle between the nth and (n + l)th amino acids 
is written as h n +\ — h n = iisin(0„i) and (p n+ 1 — ip n = 
2arcsin^i cos(0 n i)/(2i?)]. Figures [TJc) and[7Kd) investi¬ 
gate the effect of the space angle disorder on the spin 
polarization P s and (P s ) of the peptide in the presence 
of the gate voltage. It can be seen that the peak of P s -E 
is shifted towards the lower energy with the increase of D 
and the peak value nearly remains unchanged. Besides, 
(P s ) decreases quite slowly with increasing D within a 
wide range of the space angle disorder and remains large 
in the case of large D and V g . From the above, we can 
conclude that the spin filtration efficiency of the a-helical 
peptide is robust against both the on-site energy disorder 
and the space angle disorder under the gate voltage. 

IV. CONCLUSIONS 

In summary, we investigate the influence of the gate 
voltage on the quantum spin transport along the a-helical 
peptide/protein molecule contacted by two nonmagnetic 
electrodes. The spin filtration efficiency of the a-helical 
peptide/protein can be enhanced significantly by modu¬ 
lating the magnitude as well as the direction of the gate 
field, and is robust against the dephasing, the on-site en¬ 
ergy disorder and the space angle disorder in the presence 
of the gate voltage. A constructive scheme is provided for 
further experimental studies on protein spintronics and 
meanwhile, could be readily carried out and checked. 
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